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ABSTRACT 

Recently there have been claims on model-independent evidence of dynamical dark energy. Herein 
we consider a fairly general class of cosmological models with a time-evolving cosmological term of the 
form A (H) = Co + ChH 2 +C^H, where H is the Hubble rate. These models are well motivated from 
the theoretical point of view since they can be related to the general form of the effective action of 
quantum field theory in curved spacetime. Consistency with matter conservation can be achieved by 
letting the Newtonian coupling G change very slowly with the expansion. We solve these dynamical 
vacuum models and fit them to the wealth of expansion history and linear structure formation data. 
The results of our analysis indicate a significantly better agreement as compared to the concordance 
ACDM model, thus supporting the possibility of a dynamical cosmic vacuum. 

Subject headings: cosmology: observations, theory (dark energy, large-scale structure of universe) - 
methods: numerical, statistical 


1. INTRODUCTION 

The positive evidence that our Universe is speeding up 
owing to some form of dark energy (DE) pervading all 
corners of the interstellar space seems to be nowadays be¬ 
yond doubt after the first measurements of distant super¬ 
novae (Riess et al. 1998; Perlmutter et al. 1999) and the 
most recent analysis of the precision cosmological data 
by the Planck collaboration (Planck Collaboration XIII 
2015). The ultimate origin of such positive acceleration is 
unknown but the simplest possibility would be the pres¬ 
ence of a tiny and positive cosmological constant (CC) 
in Einstein’s field equations, A > 0. This framework, the 
so-called concordance or ACDM model, seems to describe 
quite well the observations (Planck Collaboration XIII 
2015) but, unfortunately, there is little theoretical mo¬ 
tivation for it. The CC is usually associated to the en¬ 
ergy density carried by the vacuum, through the param¬ 
eter pa = A/(87tG) (in which G is the Newtonian cou¬ 
pling), although it is difficult to reconcile its measured 
value (pa ~ 10 -47 GeV 4 ) with typical expectations in 
quantum field theory (QFT) and string theory, which 
are many orders of magnitude bigger. Such situation 
has triggered in the past - only to find it reinforced at 
present - the old CC problem and the cosmic coincidence 
problem (Weinberg 1989; Sahni & Starobinsky 2000; 
Padmanabhan 2003; Peebles & Ratra 2003). Both 
problems lie at the forefront of fundamental physics. 

Different theoretical scenarios have been proposed 
since long. In this Letter we take seriously the recent 
observational hint that the DE could be dynamical as 
a means to alleviate some tensions recently observed 
with the ACDM (Sahni, Shafieloo & Starobinsky 2014). 
Specifically, we focus on the dynamical vacuum mod- 

sola@ecm.ub.edu 
adriagova@ecm.ub.edu 
j decrupe7@alumnes. ub.edu 

1 High Energy Physics Group, Dept. ECM, Univ. de 
Barcelona, Av. Diagonal 647, E-08028 Barcelona, Catalonia, 
Spain 

2 Institut de Ciencies del Cosmos (ICC), Univ. de Barcelona, 
Av. Diagonal 647, E-08028 Barcelona, Catalonia, Spain 


els of the form A (H) = Co + ChH 2 + C^H, in which 
H = a/a and H = dH/dt are the Hubble rate and its 
cosmic time derivative, with Co ^ 0 a constant. We as¬ 
sume that at least one of the coefficients Ch and C^ is 
nonvanishing. Such models possess a well-defined ACDM 
limit (Ch, C^ —> 0) and involve two time derivatives of 
the scale factor, thereby can be consistent with the gen¬ 
eral covariance of the effective action of QFT in curved 
spacetime. While the general structure of A (H) can 
be conceived as an educated phenomenological ansatz, 
it can actually be related to the quantum effects on 
the effective vacuum action due to the expanding back¬ 
ground, in which the leading effects may generically be 
captured from a renormalization group equation (Sola 
2008; Shapiro & Sola 2009; Sola 2013; Sola & Gomez 
2015; Sola 2015). The dimensionless coefficients Ch and 
Ch are actually related to the /?-function of the running 
and are therefore naturally small. In the presence of 
matter conservation this is possible by letting G = G(H) 
be dynamical as well (Sola 2008, 2013). A generaliza¬ 
tion of pa(H) with higher powers of the Hubble rate, i.e. 
H n (n > 2), has been recently used to describe inflation, 
see e.g. (Lima, Basilakos & Sola 2013; Sola 2015). 

In the following we solve these dynamical vacuum mod¬ 
els A (H) and test them in the light of the recent obser¬ 
vational data, and compare their performance with the 
concordance ACDM model. 

2. BACKGROUND COSMOLOGICAL SOLUTIONS 

The field equations for the dynamical vacuum energy 
density in the Friedmann-Lemaitre-Robertson-Walker 
(FLRW) metric in flat space are derived in the standard 
way and are formally similar to the ones with strictly 
constant G and A terms: 

3H 2 =87rG(H)(p m + p r + p A (H)) ( 1 ) 

3 H 2 + 2 H= -8tt G(H) (p A + p r ) , (2) 

where pa(H) = A(H) / (8nG(H)) is the dynamical vac¬ 
uum energy density, pa(H) = -~pa(H), and G(H ) is the 
dynamical gravitational coupling. It is convenient to take 
into account from the beginning the effect of relativistic 
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TABLE 1 

Best-fit values for G 1-type models 


Model l^7l (BBN,CMB), Omh 2 (all data )v _ v _ erg ag _ y/ 2 /dof _ y^/do f _AIC AIC 


ACDM 

Yes 

0-2787S 

0.276 ± 0.004 

- 

- 

0.815 

0.815 

828.84/1010 

828.69/1010 

830.84 

830.69 

G1 

(10%,5%), Yes 

0.278 ± 0.006 

0.275 ± 0.004 

0.0015t£ml5 

0 0091 + U UUi4 

u.uuzi_ 0 0016 

0.797 

0.784 

822.82/1009 

821.97/1009 

826.82 

825.97 

ACDM 

No 

0.292 ± 0.008 

0.286 ±0.007 

- 

- 

0.815 

0.815 

583.38/604 

582.74/604 

585.38 

584.74 

G1 

(10%,5%), No 

0.290 ±0.011 

0.281 ± 0.005 

O.OOOS^.Sols 

0.0015 ±0.0014 

0.795 

0.771 

577.62/603 

575.70/603 

581.62 

579.70 

ACDM* 

Yes* 

0.297 ± 0.006 

0.293 ±0.006 

- 

- 

0.815 

0.815 

806.68/982 

806.17/982 

808.68 

808.17 

Gl* 

(10%,5%), Yes* 

0.296 ± 0.009 

0.287 ± 0.004 

0.0006 ± 0.0015 

“onoIF™ - 

0.803 

0.770 

802.66/981 

799.15/981 

806.66 

803.15 


Note. — The best-fitting values for the Gl-type models and their statistical significance (x 2 -test and Akaike information criterion AIC, see the 
text). All quantities with a bar involve a fit to the total input data, i.e. the expansion history (Omh 2 ±BAO±SNIa), CMB shift parameter, the 
indicated constraints on the value of AG/Go at BBN and at recombination, as well as the linear growth data. Those without bar correspond to a 
fit in which we use all data but exclude the growth data points from the fitting procedure. “Yes” or “No” indicates if Omh 2 enters or not the fit. 
The starred scenarios correspond to removing the high redshift point 2 : = 2.34 from Omh 2 (see text). The quoted number of degrees of freedom 
( dof ) is equal to the number of data points minus the number of independent fitting parameters. The fitting parameter v includes all data. 


TABLE 2 

Best-fit values for G2-type models 


Model It^tI (CMB), Omh 2 _ 0^ (all data) _ is e fi __ z/ e ff _ as as X 2 /d°f X 2 /d°f AIC AIC 


ACDM 

-, Yes 

0.278+S 

0.276 ±0.004 

- 

- 

0.815 

0.815 

828.84/1009 

828.69/1009 

830.84 

830.69 

G2 

5%, Yes 

0.278 ± 0.006 

0.277 ± 0.004 

0.0038iX23 

TooITSS - 

0.774 

0.773 

817.17/1008 

817.26/1008 

821.17 

821.26 

ACDM 

-, No 

0.292 ± 0.008 

0.286 ±0.007 

- 

- 

0.815 

0.815 

583.38/603 

582.74/603 

585.38 

584.74 

G2 

5%, No 

0.287 ±0.011 

0.283 ±0.005 

0.0025+"™ 

0.0030^™ 

0.763 

0.767 

572.68/602 

572.99/602 

576.68 

576.99 

ACDM* 

-, Yes* 

0.297 ± 0.006 

0.293 ± 0.006 

- 

- 

0.815 

0.815 

806.68/981 

806.17/981 

808.68 

808.17 

G2* 

5%, Yes* 

0.295 ± 0.009 

0.289 ±0.005 

0.0015+^25 

o.oo28+s:^; 

0.789 

0.765 

798.85/980 

797.05/980 

802.85 

801.05 


Note. — As in Table 1, but for G2 models with = 1 so as to maximally preserve the BBN bound (see text). The effective G2-model fitting 
parameter in this case is i/ e ff = is/4. The constraint on | AG/Go| from CMB anisotropies at recombination is explicitly indicated. 


matter, i.e. p r = (1/3 )p r , together with dust (p m = 0). 
We consider the following two realizations of the dynam¬ 
ical vacuum model: 

G1 : A(7J) = 3(co + vH 2 ) (3) 

G2: A(H,H) = 3(c 0 + vH 2 + ^aH), (4) 

o 

where we have redefined Go = 3co, Ch = 3zz and Cpj = 
2a for convenience. Model G1 is of course a particular 
case of Model G2, but it will be useful to distinguish 
between them. We can combine (1) and (2) to obtain 
the equation of local covariant conservation of the energy, 
i.e. V M (GT/ 0 ) = 0. Explicitly, since we assume matter 
conservation (meaning p m + 3 Hp m = 0 and p r + 4 Hp r = 
0), it leads to a dynamical interplay between the vacuum 
and the Newtonian coupling: 


G(pm + Pr + pa) + Gp A — 0. 


( 5 ) 


Trading the cosmic time for the scale factor a, the previ¬ 
ous equations amount to determine G as a function of a. 
Using the matter conservation equations, we arrive at 


G(a) = -Go 


« (E 2 (a)Y 

3f l m a~ 3 + 4f l r a -4 


( 6 ) 


where Go = G(a = 1) is the present value of G, 
and E(a) = H(a)/Ho is the normalized Hubble rate 
to its present value. The prime stands for d/da , and 
fli = Pm/PcO (with p c o = 3 Hq /8ttGq) are the currently 
normalized energy densities with respect to the critical 
density. Inserting (4) and the above result for G(a) in 


Eq. (1) and integrating, we obtain: 


E 2 {a) = 1 + 


- - i; 


-1 + «- 4 «' 



( 7 ) 


where have introduced 

C = = 1 - Gff ; ^ = 2—= 1 - • (8) 

1 — a 1 — 


For small \v,a\ 1 (the expected situation), we can use 

the approximations v e Q ~ v — a and i/ ff ~ v — (4/3)a. 
Note that, in order to simplify the presentation, we have 
removed from Eq. (7) terms proportional to H r -C 
that are not relevant here. We can check e.g. that 
in the radiation-dominated epoch the leading term in 
the expression (7) is ~ fl r a~ 4 ^ , whilst in the matter- 
dominated epoch is ~ H m a -3 ^. Furthermore, we find 
that the (full) expression for E 2 (a) reduces to the ACDM 
form, 1 + H m (a -3 — 1) + H r (a -4 — 1), in the limit i/, a —> 0 
(i.e. £,£' —>• 1). Notice also the constraint among 

the parameters, Co = Hq [Ry — v + a (f2 m + | H r )] , 
which follows from matching the vacuum energy density 
Pa{H) to its present value p 3 for H = Hq and using 
+ n r + Ha = 1. The explicit scale factor dependence 
of the Newtonian coupling ensues upon inserting (7) in 
(6) and computing the derivative. We refrain once more 
from writing out the full expression here, but one can 
check that in the limit a —» 0 (relevant for the Big Bang 
Nucleosynthesis epoch) it behaves as 


G(a) = G 0 a 4(1 - ?,) cs G 0 (l + 4i/' ff In a). (9) 

Thus, the gravitational coupling evolves logarithmically 
with the scale factor and hence changes very slowly. This 
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FlG. 1.— Left: Evolution of the effective EoS ojd(z ), Eq. (11), for the models under consideration: Right: The corresponding evolution of the 
effective DE density Qd(z) normalized to the critical density (see text). 



logarithmic law was motivated previously in (Sola 2008, 
2013) within the context of the renormalization group 
of QFT in curved spacetime. For is = a = 0 we ob¬ 
tain G = Go identically, i.e. the current value of the 
gravitational coupling. However the situation G = Go 
is also attained in the limit a —0 for v = (4/3)a (i.e. 

= 1); and indeed we shall adopt this setting here¬ 
after in order to maximally preserve the BBN constraint 
for the G2 model. The effective fitting parameter will 
be is e f{ = v/A. Obviously this setting is impossible for 
Gl, so in this case we will adopt the average BBN re¬ 
striction |AG/G| < 10% in the literature (Chiba 2011; 
Uzan 2011). At the same time we require |AG/G| < 5% 
at recombination (z ~ 1100) for both Gl and G2 from 
the CMB anisotropy spectrum (Chiba 2011). 

The expression for the dynamical vacuum energy den¬ 
sity can be obtained from Friedmann’s equation (1), in 
combination with the explicit form of G(a). We quote 
here only the simplified expression valid for the matter- 
dominated epoch: 


Pa(«) — PcO a 



+ ^(1-{-<>«) 


( 10 ) 


For £ — >• 1 we have pa —t p c o(l — D m ) = PcO^A and we 
retrieve the ACDM case with strictly constant pa- The 
form (10) is sufficient to obtain an effective DE density 
Pd(z) an d effective equation of state (EoS) for the DE at 
fixed G = Go, as conventionally used in different places 
of the literature - see e.g. (Sola & Stefancic 2006, 2005; 
Shafieloo et al. 2006; Basilakos & Sola 2014). We find: 

wj)(a) =- , , „ . (11) 

i I pm\Q -) G(fl) Go 

1 + PA(a) G(a) 

Iii Fig. 1 (left) we plot ujd as a function of the cosmic 
redshift z = — 1 + 1/a for models Gl and G2. Near our 
time, usd stays very close to —1 (compatible with the 
ACDM), but at high z it departs. In the same Fig. 1 
(right) we plot Qd{z) = Pd{z)/ p c {z), i.e. the normal¬ 
ized DE density with respect to the critical density at 
constant Go- The asymptotes of wb for each model at 
z > 4 are due to the vanishing of f!o(z) at the corre¬ 
sponding point (as clearly seen in the figure)- confer the 
aforementioned references for similar features. 


3. FITTING THE MODELS TO THE OBSERVATIONAL 
DATA 

Let us now test these models versus observation. First 
of all, we use the available measurements of the Hub¬ 
ble function as collected in (Ding et al. 2015). These 
are essentially the data points of (Farooq & Ratra 2013) 


in the redshift range 0 St z ^ 1.75 and the BAO mea¬ 
surement at the largest redshift H(z = 2.34) taken after 
(Delubac 2015) on the basis of BAO’s in the Lyct for¬ 
est of BOSS DR11 quasars. We define the following \ 2 
function, to be minimized: 


Xomh 2 — 


= E E 


i= 1 j=i -\-1 


Omht h (Hj, Hi) - Omh 2 obs (Hi, Hj) 


®Omh 2 i,j 


( 12 ) 

where N is the number of points H(z) contained in 
the data set, Hi = H(z-i), and the two-point diagnostic 
Omh 2 (z 2 , zi) = [h 2 (z 2 ) - /i 2 (,zi)]/[(l + z 2 ) 3 - (1 + ^i) 3 ] 
was defined in (Sahni, Shafieloo & Starobinsky 2014), 
with h(z) = hE(z), and <Jo m h 2 i,j is the uncertainty 
associated to the observed value Omhl bs (Hi, Hj) for a 
given pair of points ij, viz. 


®Omh 2 i.j 


4 


h2 ( Z i) a h( Zi ) + h2 ( Z o) a h( Zj ) 

[{l + Zi) 3 ~ (1 + Zj) 3 } 2 


(13) 


For the ACDM the two-point diagnostic boils down to 
Omh 2 (z 2 , z i) = h 2 , which is constant for any pair z i, 
Z 2 - Using this testing tool and the known observational 
information on H{z) at the three redshift values 2 = 
0,0.57, 2.34 the aforementioned authors observed that 
the average result is: Omh 2 = 0.122 ± 0.010, with very 
little variation from any pair of points taken. The ob¬ 
tained result is significantly smaller than the correspond¬ 
ing Planck value of the two-point diagnostic, which is 
constant and given by Omh 2 = h 2 = 0.1415 ± 0.0019 
(Planck Collaboration XIII 2015). 

A departure of Omh 2 from the Planck result should, 
according to (Sahni, Shafieloo & Starobinsky 2014), sig¬ 
nal that the DE cannot be described by a rigid cosmo¬ 
logical constant. For the ACDM we obtain Omh 2 = 
0.1250 ± 0.0039, and Omh 2 = 0.1402 ± 0.0059, by taking 
all data points, and excluding the high redshift one, re¬ 
spectively. Since there is a priori no reason to exclude the 
high-redshift point (Delubac 2015), whose uncertainty is 
one of the lowest in the full data sample, relaxing the ten¬ 
sion with data may require the dynamical nature of the 
DE. The vacuum models Gl and G2 considered here, 
Eqs. (3,4), aim at cooperating in this task. 

For these models the theoretical value Omh 2 h of the 
two-point diagnostic entering (12) can be computed, in 
the matter-dominated epoch (relevant for such observ¬ 
able), as follows: 


Omh 2 G (zi, Zj) 


h 2 (1 + Zi ) 3 ^ — (1 + XJ )« 
~T~ (1 + Zi) 3 — (1 + Zj) 3 


(14) 






















4 



FlG. 2.— Likelihood contours in the (Sl m , z'eff) plane (for the values 
—2 In C/Cmax — 2.30, 6.16,11.81, corresponding to la, 2a and 3a 
confidence levels for the G2 model using the full data analysis indicated 
in Table 2. The i/ e ff — 0 region (ACDM) is disfavored at ~ 3a. 

It is evident that for £ = 1 we recover the ACDM 
result, which remains anchored at Omh 2 (zi,Zj) = 
fl m h 2 ( \/Zi,\/zj ). However, when we allow some small 
vacuum dynamics (meaning v and/or a different from 
zero) we obtain a small departure of £ from 1 and there¬ 
fore the DE diagnostic Omh 2 deviates from Q m h 2 . In 
this case Omh 2 evolves with cosmic time (or redshift). 

To the above Hubble parameter data we add the re¬ 
cent supernovae type la data (SNIa), the Cosmic Mi¬ 
crowave Background (CMB) shift parameter, the Bary- 
onic Acoustic Oscillations (BAO’s), the growth rate for 
structure formation (see next section) and the BBN and 
CMB anisotropy bounds. Contour lines for v e R = 1 — £ 
are shown in Fig. 2 for model G2 at fixed /' = 1. 
The x 2 functions associated to SNIa distance modulus 
p(z), the BAO A-parameter and the CMB shift param¬ 
eter can be found in (Gomez-Valent, Sola & Basilakos 
2015). Therein, one can also find the corresponding ref¬ 
erences of the data sets that we have used in the present 
analysis. 


a (when non-relativistic matter dominates over the vac¬ 
uum) we have 5 m (a) = a s , where s = 3£ — 2 = 1 — 3zz e g. 
If £ = 1 [v e s = 0), then 5 m (a ) ~ a and we recover 
the ACDM behavior. Thus, the initial conditions set at 
a high redshift = (1 — ai)/aj, say = 100 (or at 
any higher value), are the following. For the growth 
factor we have 5 m (ai) = af, and for its derivatives: 

6m( a i) = sa i~\ C(ai) = s(s - 1K" 2 . 

In practice we investigate the agreement with the 
structure formation data by comparing the theoretical 
linear growth prediction f(z) = —(1 + z)d\nS m /dz and 
the growth rate index 'y(z) with the available growth 
data - following (Gomez-Valent, Sola & Basilakos 2015) 
and references therein. Recall that 7 is defined 
through f{z) ~ f 2 m (z) 7 ( 2 ), and one typically ex¬ 
pects 7 ( 6 ) = 0.56 ± 0.05 for ACDM-like mod¬ 

els (Pouri, Basilakos & Plionis 2014). A most con¬ 
venient related quantity is the weighted growth rate 
f(z)a$(z), cf. (Song & Percival 2009), where ag{z) is the 
rms mass fluctuation amplitude on scales of Rg = 8 h ” 1 
Mpc at redshift z. The latter is computed from 

S m (z) \ f 0 °° k n ‘+ 2 T 2 (n m ,k)W 2 (kR 8 )dk l 1 
^ CT8 ’ A ^(0) [ k n ’>+ 2 T 2 (n rntA ,k)W' 2 (kR 8 )dk 

(16) 

W being a top-hat smoothing function (see e.g. 
(Gomez-Valent, Sola & Basilakos 2015) for details) 
and T(S2 m ,fc) the transfer function, which we take 
from (Bardeen et al. 1986). The values of cr$ = CTg(0) 
for the various models are collected in Table 1, and in 
Fig. 3 we plot f(z)cr$(z) and 7 (z) for them. 

The joint likelihood analysis is performed on the set 
of Omh 2 +BAO-|-SNIa-|-CMB, BBN and linear growth 
data, involving one (f 2 m ) or two (fl m , v e g) independently 
adjusted parameters depending on the model. For the 
ACDM we have one parameter ( n p = 1) and for G1 and 
G2 we have n p = 2. Recall that for G2 we have fixed 

e = i- 


4. LINEAR STRUCTURE FORMATION 

Finally, we take into consideration the data on the lin¬ 
ear structure formation. For the G1 and G2 models the 
calculation of S m = Sp m /p m is significantly more compli¬ 
cated than in the ACDM case and follows from applying 
linear perturbation theory to Einstein’s field equations 
and Bianchi identity (5) (Grande et al. 2010). The final 
result reads 3 : 

X" qA' 

C+^Q-^m)+^(8-nn m +3n 2 m -aQ' m ) = o, 

(15) 

with flm(a) = 87t G(a)p m (a)/3H 2 (a). Notice that the 
(iz,a) model-dependence is encoded in H(a) - cf. Eqs. 
(7,8). To solve the above equation (numerically) we have 
to fix the initial conditions for S m , S' m and <5" . We take 
due account of the fact that for these models at small 

3 The third-order feature of this equation is characteristic of the 
coupled systems of matter and DE perturbations for cosmologies 
with matter conservation, after eliminating the perturbations in 
the DE in favor of a single higher order equation for the matter 
part — cf. (Gomez-Valent, Karimkhani & Sola 2015) for details. 
For A =const. Eq. (15) boils down to the (derivative of the) second 
order one of the ACDM. 


5. DISCUSSION 

The main results of this work are synthesized in Ta¬ 
bles 1-2 and Figures 1-3. In particular, from Fig. 2 
we see that the model parameter is e ff for G2 is clearly 
projected onto the positive region, which encompasses 
most of the 3er range. Remarkably, the x 2 - v alue of the 
overall fit is smaller than that of ACDM for both G1 
and G2 (cf. Tables 1-2). To better assess the distinc¬ 
tive quality of the fits we apply the well known Akaike 
Information Criterion (AIC) (Akaike 1974), which re¬ 
quires the condition Nt a t/n p > 40 (amply satisfied in 
our case). It is defined, for Gaussian errors, as follows: 
AIC = —2 In £ max + 2n p = xLn + 2 n p> where £ max (resp. 
Xmi n ) is the maximum (resp. minimum) of the likeli¬ 
hood (resp. x 2 ) function. To test the effectiveness of 
models Mi and Mj , one considers the pairwise difference 
(AAIC)y = (AIC)j — (AIC)j. The larger the value of 
Ay = |A(AIC)y|, the higher the evidence against the 
model with larger value of AIC, with Ay > 2 indicating 
a positive such evidence and Ay > 6 denoting significant 
such evidence. 

From Tables 1-2 we see that when we compare the fit 
quality of models i =G1, G2 with that of j = ACDM, 
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FlG. 3.— Left: Comparison of the observed data with error bars (in green) and the theoretical evolution of the weighted growth rate of clustering 
f(z)as(z) for each dynamical vacuum model and the ACDM. Right: The corresponding evolution of the linear growth index 7 (z). 


in a situation when we take all the data for the fit opti¬ 
mization, we find A y - = AAIC > 9.43 for G2 and 4.72 
for Gl, suggesting significant evidence in favor of these 
models (especially G2) against the ACDM - the evidence 

ratio (Akaike 1974) being ER = e Si ’/ 2 > 111.6 for G2 
and 10.6 for Gl. Worth noticing is also the result of 
the fit when we exclude the growth data from the fitting 
procedure but still add their contribution to the total % 2 . 
This fit is of course less optimized, but allows us to risk 
a prediction for the linear growth and hence to test the 
level of agreement with these data points (cf. Fig. 3). 
It turns out that the corresponding AIC pairwise differ¬ 
ence with the ACDM are similar as before (cf. Tables 
1-2). Therefore, the ACDM appears significantly disfa¬ 
vored versus the dynamical vacuum models, especially in 
front of G2, according to the Akaike Information Crite¬ 
rion. Let us mention that if we remove all of the H{z) 
data points from our analysis the fit quality weakens, 
but it still gives a better fit than the ACDM (cf. the 
third and fourth row of Tables 1 and 2). If, however, 
we keep these data points but remove only the high red- 
shift point z = 2.34 (Delubac 2015), the outcome is not 
dramatically different from the previous situation (con¬ 
fer the starred scenarios in Tables 1 and 2), as in both 
cases the significance of v b q ^ 0 is still close to ~ 2 a 
with A ij > 7 for G2 (hence still strongly favored, with 
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